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Abstract — This paper investigates the propagation characteris- 
tics of circular waveguides whose interior surface is coated with 
a thin metamaterial liner possessing dispersive, negative, and 
near-zero permittivity. A field analysis of this system produces 
the dispersion of complex modes, and reveals in detail intriguing 
phenomena such as backward-wave propagation below the un- 
lined waveguide's fundamental-mode cutoff, resonant tunneling 
of power, field collimation, and miniaturization. It is shown how 
the waveguide geometry and metamaterial parameters may be 
selected to engineer the lined waveguide's spectral response. The- 
oretical dispersion and transmission results are closely validated 
by full-wave simulations. 

Index Terms — Metamaterials, circular waveguides, inhomo- 
geneous waveguides, epsilon-near-zero, negative permittivity, 
backward-wave, below-cutoff propagation, miniaturization 

I. Introduction 

HOLLOW waveguides are widely used in applications 
requiring high power-handling capability and simple 
integration with radiating devices such as horn antennas. 
As such, they are ideally suited to microwave radio links, 
radar, and satellite communications. It is well known that 
air-filled waveguides enclosed by perfect-electric-conducting 
(PEC) boundaries support a discrete spectrum of modes, each 
possessing a cutoff frequency, f c , above which it is capable of 
propagating power. Homogeneously filling the vacuum region 
of the waveguide with an isotropic dielectric material serves 
to reduce these cutoff frequencies, without modifying the 
corresponding transverse modal field distributions. However, 
it is sometimes desirable to have access to the enclosed region 
of a waveguide, making it convenient to partially fill it with 
a dielectric to achieve similar results. These inhomogeneous 
structures have special impedance-boundary conditions and 
support modal field distributions that can differ significantly 
from those of their conventional, homogeneously-filled coun- 
terparts. 

One such instance is the PEC circular waveguide containing 
two concentric dielectric regions. This structure has been 
extensively investigated in the form of a waveguide loaded 
with a dielectric rod (TJ, (2). However, the same analysis can 
be applied to its variant: a waveguide whose interior wall 
is lined with a dielectric medium. This structure supports 
hybrid electric (HE) and hybrid magnetic (EH) modes which 
are similar to the transverse electric (TE) and transverse 
magnetic (TM) modes of a homogeneously filled circular 
waveguide, except that the longitudinal electric and magnetic 
fields do not, in general, vanish. In certain cases, such as 
circularly symmetric modes or near cutoff, the dielectric-lined 
circular waveguide's HE and EH modes can be likened to 



the TE and TM modes, respectively, of its homogeneously- 
filled counterpart. The inhomogeneous nature of the waveguide 
volume introduces additional degrees of freedom in engineer- 
ing its modal dispersion characteristics, such as the dielectric 
layer's permittivity, permeability, and thickness, which may 
be used to enhance bandwidth and introduce backward-wave 
propagation (2), (5). 

Metamaterial loading of waveguides is a new avenue of re- 
search in which the propagation characteristics of an unloaded 
waveguide are altered through the inclusion of metamaterials, 
whose characteristics may be engineered to achieve intriguing 
propagation phenomena. For example, the theory of images 
enables propagation inside an infinite array of metamaterial 
inclusions to be described by placing a finite number of such 
inclusions inside a waveguide. In fact, it has been shown that 
propagation can be restored in below-cutoff waveguides by 
loading them with chains of magnetic and/or electric scatterers, 
such as the split-ring resonator (SRR), complimentary split- 
ring-resonator (CSRR), or wire-lines (4), (5). In these works, 
the metamaterial-loaded waveguides are often treated as ho- 
mogeneously filled by a material possessing epsilon-positive 
(EPS), mu-positive (MPS), epsilon-negative (ENG), and/or 
mu-negative (MNG) responses. Simultaneous EPS and MPS 
is referred to as 'double-positive' (DPS), and simultaneous 
ENG and MNG is referred to as 'double-negative' (DNG). 
However, whereas homogeneously filled models may explain 
certain elementary dispersion features, they cannot describe 
the full dispersion of complex modes, the distribution of 
power in the transverse cross-section, and the possibility 
of backward coupling of power, all of which have been 
observed in inhomogeneously filled structures. Furthermore, 
these homogeneous models prove inadequate in applications 
in which the waveguide interior must remain empty. Therefore, 
a more complete description of propagation phenomena in 
metamaterial-loaded waveguides demands an inhomogeneous 
model. 

In 2004, Alu et al observed that the propagation of TE 
and TM modes can be made independent of the thickness 
of a parallel-plate waveguide, provided that it is filled inho- 
mogeneously with bilayers of the varieties DPS -DNG, DNG- 
ENG, DPS-MNG and ENG-MNG (6). It was later shown by 
Silveirinha and Engheta that complete transmission or 'su- 
percoupling' could be achieved in discontinuous waveguides 
of arbitrarily small cross section, provided that they are ho- 
mogeneously filled with single-negative (SNG) metamaterials 
possessing epsilon-near-zero (ENZ) values (T). Subsequent 
experimental validations of supercoupling employed either a 
waveguide containing an embedded metamaterial with zero 
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Fig. 1. Transverse cross-section of the dielectric -lined circular waveguide 
consisting of two concentric dielectric regions. An inner core region of radius 
a and a liner of thickness b — a are filled with materials described by relative 
parameters {e r i, /iri} and {e r 2, Mr2}> respectively. 



permittivity, such as the CSRR [8], or a narrow waveguide 
channel that naturally exhibits zero effective permittivity at its 
cutoff frequency (9). 

Whereas the standard analysis of the dielectric-lined circular 
waveguide has been applied to the case of metamaterial 
liners (TO), (TTj, these works have either omitted the necessary 
frequency dispersion, concerned themselves with the DNG 
case alone, or limited their study to obtaining the field expres- 
sions and dispersion relations. The present work investigates 
the intriguing effects and implications of partially filling 
circular waveguides with thin, dispersive ENZ metamaterial 
liners and establishes several phenomena, including below- 
cutoff propagation, resonant tunneling and backward coupling 
of power, field collimation, and the possibility of waveguide 
miniaturization. Section [TT] presents a rigorous hybrid-mode 
analysis to determine the lined waveguide's general dispersion 
features, cutoff frequencies, and field configurations. Sec- 
tion III presents a representative metamaterial-lined waveguide 
analyzes its dispersion features, 



design, analyzes its dispersion features, and offers design 
guidelines for engineering the frequency response of its HEn 
mode. Finally, Sec. |IV] presents full-wave dispersion and 
transmission simulations showing close agreement with the 
theoretical results. Potential applications of this study include 
the design of metamaterial-lined waveguides to demonstrate 
reverse Cherenkov radiation, as in Ref. (12), and the inclu- 
sion of metamaterials inside magnetic-resonance scanners to 
enable traveling- wave imaging at low static field strengths, as 
proposed in Ref. |T3) . 

II. Theory 

Figure [T] presents the geometry of the dielectric-lined cir- 
cular waveguide under consideration. An inner core region 
of permittivity e\ = e r ieo and permeability \i\ = p r iPo is 
surrounded by a dielectric material of thickness t — b — a with 
permittivity e 2 = e r 26o and permeability p 2 = Mr 2 Mo- The 
dielectric-lined waveguide modes are no longer pure TE mn and 
TM mn modes; they are, instead, classified as HE mn and EH mn 
modes generally possessing both H z ^ and E z ^ 0. We now 
present the analysis of the dielectric-lined circular waveguide, 
which may also be found in several prior works (Tj> j2j, (TO), 
in order to provide the reader with some context for the main 
analytical results of this work. Suppressing the propagation 
term (e~ lz ) and the time-harmonic nature of the fields (e Ju;£ ), 



the electric- and magnetic-field components in the axial (z-) 
direction in region 1 (0 < p < a) may be represented by 



Eiz = CiJ n (ki p p) cos(n</>) 
H lz = C 2 J n {k lp p) sin(ra0) 



(la) 
(lb) 



where C\ and C 2 are the (generally complex) amplitude 
coefficients, J n (ki p p) is a Bessel function of the first kind, n 
is the azimuthal mode index, k\ p = k\ — 7 2 , k\ = w(eipi) - 5 , 
and 7 is the axial propagation constant. In region 2, the 
radially outgoing and incoming fields are represented by 
Bessel and Neumann functions, J n (k2 P p) and Y n (k2 P p). Since 
the tangential electric-field components must vanish at the 
PEC boundary (p — b), the solution in region 2 (a < p < b) 
can be shown to have the following form: 



E 2z = C 3 F n (k 2p p) cos(n0) 
H 2z = C A G n (k 2p p) sin(n0) 



(2a) 
(2b) 



where 



F n (k 2p p) = Y n (k 2p b)J n (k 2p p) - J n {k 2p b)Y n (k 2p p) (3a) 

G n (k 2p p) = Y;(k 2p b)J n (k 2p p) - J' n (k 2p b)Y n (k 2p p) (3b) 

Here, C3 and C4 are the amplitude coefficients, k\ p = k\—rf, 
and k 2 = w(e 2 p 2 ) - 5 . Knowing the longitudinal field compo- 
nents, the tangential field components can be obtained using 
the transverse decomposition of Maxwell's equations. The 
transverse components of the electric and magnetic fields (that 
is, the p and (j) components) in both regions are given by 
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where 

F' n {k 2pP ) = Y n {k 2p b)J' n {k 2p p) - J n {k 2p b)Y^k 2pP ) (5a) 

G' n (k 2p p) = Y^{k 2p b)J' n {k 2pP ) - J' n {k 2p b)Y^k 2pP ) (5b) 

The continuity of the tangential electric- and magnetic- field 
components at p = a relate the coefficients to one another as 
follows: 
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where equation (6b) is used for EH modes and equation 
(6c) is used for HE modes. Using equations (6b) and (6c), 
the coefficients C\ and C 2 can be eliminated to obtain the 
following dispersion relation: 

1 1 
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where 



B = [ 



ei J' n (k lp a) _ e 2 F^(k 2p a) 
ki p J n (k lp a) k 2p F n (k 2p a) 

Mi J' n {k lp a) n 2 G' n {k 2p a) 



(7a) 



(7b) 



(7c) 



ki p J n {ki p a) k 2p G n {k 2p a) 

Dielectric-lined waveguides are inhomogeneous structures that 
support hybrid modes with complex propagation constants of 
the form 7 = a + jf3. Whereas nonzero a is attributed to 
dielectric or conductive loss in homogeneous waveguides, its 
existence in inhomogeneous waveguides is a direct result of 
the coupling of power between regions and can appear even 
in the absence of losses. Modes for which the value of a 
is substantially larger than f3 are termed evanescent. These 
modes attenuate rapidly in the longitudinal direction, and are 
limited in their ability to transport power along the waveguide, 
much like the evanescent modes of a homogeneously filled 
waveguide; as a result, they shall not be the focus of this 
study. On the other hand, modes with a much smaller than 
/3 can lead to the propagation of real power; we term these 
propagating modes. Modes with a comparable in magnitude to 
/3 can be considered neither purely propagating nor evanescent, 
and are therefore termed complex modes. These complex 
modes exist in complex conjugate pairs, and the two conjugate 
modes transport power in opposite directions. Only if the 
conjugate modes are excited unequally can they transport real 
power |T4| . 

Using Muller's method to determine the roots in the 
complex-7 plane, it is possible to numerically arrive at a set 
of continuous dispersion curves. It should be noted that there 
exist certain limiting cases in which the modal fields revert to 
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Fig. 2. The HE\ \ -mode cutoff frequency versus liner permittivity. The full 
dispersion relations of equation (7a) (solid black line) is compared to the 
approximate expression equation (8a) (empty grey circles). The waveguide's 
dimensions are b = 15mm and a = 14mm (t = 1mm). The permittivity of 
the liner e r 2(w) (dashed grey line) follows a Drude model with parameters 
w ep = 3.550GHz and uot = 0. 



those of the homogeneous case. For 7 = 0, (i.e., at cutoff) the 
RHS of the dispersion relation (7a) becomes zero and the HE 
and EH modes are decoupled. This results in the roots of (7b) 
and (7c), determining, respectively, the cutoff frequencies of 
the EH and HE modes. Since the fundamental TEn mode is 
the most widely used in applications, this paper focuses on 
the propagation characteristics of the HEu mode. 

In Fig. [2| we show the relationship between the cutoff 
frequency of the HEu mode and the liner's relative per- 
mittivity, e r2 (solid black line), which is obtained using the 
analytically derived dispersion relation (7c) through enforcing 
the cutoff condition, 7 = 0. The representative circular 
waveguide investigated has physical dimensions b = 15mm 
and a = 14mm (liner thickness of t = 1mm), with an inner 
vacuum region (e r i = /i r i = 1) and a liner region with 
dispersive e r2 and a nonmagnetic response (fi r2 = 1). It should 
be noted that e r2 is here an independent variable; i.e., each 
point on the curve is only concerned with e r2 at a single 
frequency (specifically, cutoff). Figure [2] may be employed in 
choosing the value of permittivity required to achieve a desired 
cutoff frequency, after which a suitable dispersive permittivity 
function satisfying this condition may be selected. 

A well known result of the homogeneously filled dielectric 
waveguide is that the cutoff frequency of the TEn mode 
scales as the inverse square of permittivity, when that per- 
mittivity is assumed positive and nondispersive. Even though 
the metamaterial liner only occupies a small portion of the 
total cross-section of the waveguide, it is to be expected that 
a thin EPS liner with permittivity larger than unity marginally 
lowers the natural cutoff frequency, and conversely, that an 
EPS liner with permittivity smaller than unity will increase the 
natural cutoff frequency. Indeed, Fig. [2] shows that the cutoff 
frequency is weakly dependent on large positive permittivity 
values. However, the cutoff frequency increases dramatically 
for small positive ENZ values, suggesting that the waveguide 
is thrust more deeply into cutoff as e r2 — » + . In this region, 
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the liner permittivity is characterized by both EPS and ENZ, 
which we more compactly refer to as 'epsilon positive and 
near zero' (EPNZ). 

The cutoff frequency's dependence on negative permittivity 
values is even less intuitive: large, negative liner permittivities 
in the range shown produce marginally decreased cutoff fre- 
quencies; however, as e r 2 — » 0~, the cutoff frequency of the 
HEn mode is dramatically lowered, suggesting that a waveg- 
uide lined with a metamaterial possessing such permittivity 
values may support propagation at very low frequencies, well 
below its natural cutoff. It is evident from Fig. [2] that this 
regime is characterized by both ENG and ENZ, which we term 
'epsilon negative and near zero' (ENNZ). Although not shown, 
it has been determined that thinner liners require more extreme 
ENNZ values (i.e., closer to zero) to realize the same reduction 
in cutoff frequency. As may be surmised from the curve, there 
exists an asymptote as e r 2 — » 0~ at which the cutoff frequency 
is reduced to dc. Of course, this is a pathological case, since it 
would require an ENZ material at dc and, in any event, could 
only result in the trivial electrostatic field solution. 

Knowing that ENNZ liners reduce the cutoff frequency 
of the HEu mode, the decoupled dispersion relation for 
the HE modes (7c) can be simplified using small- argument 
approximations to reduce the Bessel and Neumann functions 
with the goal of obtaining a non-transcendental expression. 
The two conditions provided by the ENNZ liner, namely 
that f c — > and e r 2 — » 0~, validate the assumption {k2 P b, 
ki p a} — » in this frequency regime. Therefore, a first- 
order power- series approximation of the Bessel and Neumann 
functions can be applied to (7c), which yields the result that 
the cutoff frequency of the HEu mode is related to the 
permittivity, e r 2, of the liner and thickness, t = b — a, as 
follows: 

2(2)°-*c l-K . 5 

/c - ( ^^r ){ i^m ] (8a) 

where 

l + £ 

K = e r2 — (8b) 

1 a 2 

This simple relationship constitutes a design equation, since 
it specifies the ENNZ values of e r 2 necessary to produce 
a desired f c , given waveguide dimensions a and b. In Fig. 
[2j the approximate expression (8a) (empty grey circles) is 
compared to the exact analytical expression (7a) (solid black 
line). It is evident that (8a) models the trend over all e r 2 values 
with good accuracy, but is expectedly most accurate in the 
ENNZ region in which the asymptotic limits are valid and the 
cutoff frequency is strongly reduced. The margin of error in 
this expression can be evaluated by applying the two cases 
describing an unloaded waveguide: one in which e r 2 — » 1 
and simultaneously a — » 0, and the other in which a — » b, 
and comparing the results to the well-known expression for 
the 77^ i -mode cutoff frequency of a homogeneously vacuum- 
filled circular waveguide (f c ,uniined = 1.841c/(27ra)). It is 
easily shown that the former case produces f c = 2c/(27ra), 
whereas the latter case yields f c = 1.633c/ {2ira) \ the geomet- 
ric average of these limiting cases is within 2% of f c ,uniined- 



Equation (8a) also seems to suggest that the cutoff frequency 
can be set to zero by choosing K = 1 which, from (8b), yields 
the following result: 

_ 1 - b 2 /a 2 

e r 2,ma X - j + ^2 ( 9 > 

Equation (9) suggests that there exists a particular maximum 
ENNZ value of e r 2 that reduces the cutoff frequency to zero, 
which is solely dependent on the waveguide's and liner's 
physical dimensions. Inserting the above representative waveg- 
uide dimensions into (9), one obtains e r 2, ma ^ = —0.06889 
which coincides exactly with the value of permittivity at which 
the asymptote occurs in Fig. [2j as expected. Thus, the liner 
needs only possess a sufficiently ENNZ e r 2 to effect a desired 
reduction in the iffin-mode cutoff frequency. For practical 
metamaterial technologies, operation in this frequency-reduced 
regime offers the potential for miniaturization of waveguide 
components. 

III. Design 

In Fig.[2| it was implicitly assumed that the liner permittivity 
required to produce a desired cutoff frequency could actually 
be achieved at that frequency. However, designing for a 
particular //En -mode cutoff requires the dispersive nature of 
e r 2 to be taken into account. Hence, the liner permittivity's 
dispersion is now described by a Drude model with e r 2{oo) = 
l—uo1 p /(jj(uo — juo t ), in which uj ep is the plasma frequency and 
cut is the damping frequency establishing the liner's loss. This 
model accurately approximates the complex, dispersive nature 
of the liner's permittivity over a select frequency range. This 
type of dispersion can be realized by any number of existing 
metamaterial technologies including wire-grid medi a |T5) or 
radially arranged transmission line metamaterials |16|, and 
therefore lends itself readily to practical implementation. To 
design the metamaterial-lined circular waveguide for a specific 
cutoff frequency f c , the waveguide dimensions (outer radius b, 
inner radius a) are chosen and the corresponding required e r 2 
may be determined from equation (8a). Thereafter, it remains 
only to choose an appropriate dispersion model for the liner 
permittivity such that e r 2 assumes the needed value at the 
desired cutoff frequency. Conversely, the dimensions of the 
waveguide can be determined if given an e r 2(uj) based on a 
known metamaterial frequency response. 

For the chosen waveguide dimensions (b = 15mm and 
a = 14mm), it can be observed from Fig. [2] that a liner 
permittivity of e r 2 = —0.09 will result in a reduced cutoff 
frequency of f c = 3.381GHz. To achieve this goal, the Drude 
model parameters for the liner are set to u ep = 3.550GHz 
and oo t = 0MHz. Superposing the metamaterial 's dispersive 
permittivity (dashed grey curve) on the cutoff-frequency curve, 
as also shown in Fig. [2j two intersections are observed: 
the lower-frequency intersection corresponds to the designed 
HEu -mode cutoff, which, for e r 2 = —0.09, indeed produces 
f c = 3. 381 GHz; the higher-frequency intersection occurs 
at e r2 = 0.6438, which results in f c = 5.958GHz - a 
slight increase as compared to the cutoff frequency of the 
homogeneously vacuum-filled waveguide. Since each of these 
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Fig. 3. (a) Dispersion of a/ko and (3/ko for the metamaterial-lined 
waveguide's HEu mode (dashed and solid black lines, respectively) as 
compared to the dispersion of the TE\ i mode of an unlined waveguide of the 
same outer radius (empty and solid grey circles, respectively). The normalized 
complex electric-field magnitude and vectors in the transverse plane obtained 
from the theory in Sec. [n] for the lined waveguide's at (b) fs = 5.958GHz 
and (c) fi = 3.381GHz, corresponding to the respective cutoff frequencies 
of the forward- wave and backward- wave bands. 



intersections corresponds to a cutoff frequency of the HEu 
mode, it is evident that the lined waveguide supports dual- 
band operation (this is verified by transmission simulations 
presented in the next section). These cutoff frequencies may 
be adjusted simply by varying the liner's permittivity profile 
and/or the waveguide's dimensions using (8a) and/or a figure 
equivalent to Fig. [2] 

Figure [3(a)] presents the iffin-mode dispersion of the loss- 
less metamaterial-lined waveguide as obtained from equations 
(7a)-(7c). Here, the solid and dashed lines, respectively, indi- 
cate /3/ko and a/ko. For clarity, and consistent with the mode 
classification discussed in the previous section, a spectrum 
of strongly evanescent modes whose cutoffs lie below the 
frequency region shown have been omitted. Also shown for 
comparison is the T^n-mode dispersion of a homogeneously 
vacuum-filled (unlined) waveguide of the same outer di- 
mension, where /3/ko and a/ko are represented using grey 
solid and open circles, respectively. As suggested from the 
two intersections in Fig. [2] the metamaterial-lined waveguide 
supports two propagating bands: a forward- wave band near the 



cutoff of the unlined waveguide, and a frequency-reduced band 
well below this cutoff frequency. The frequency -reduced band 
has two regions of interest: a propagating backward- wave band 
for f < fi = 3. 381 GHz and a complex-propagation band 
for f 1 <f<f 2 = 3.644GHz. Backward- wave propagation 
has previously been observed in inhomogeneous waveguides, 
albeit with conventional positive- valued dielectrics, and result- 
ing in only a marginal reduction in cutoff frequency (T7). 
The use of ENNZ materials in the present work suggests that 
backward-wave propagation can be achieved at dramatically 
reduced frequencies. As will be shown in the next section, the 
complex-propagation band does not allow for propagation of 
real power, a fact that may be attributed to the simultaneous 
existence of it's complex conjugate, which is excited equally 
in this region |T8| . The stopband in the region f 2 < f < fs = 
5.958GHz is explained by the fact that the metamaterial liner 
exhibits positive permittivities approaching unity as frequency 
is increased; as a result, propagation in the lined waveguide 
at these frequencies strongly resembles that in the unlined 
waveguide under its natural TEu cutoff. Accordingly, the 
forward-wave band in the region / > is also very similar 
in curvature and shape to that of the unlined waveguide. 
Furthermore, as evident from Fig. [2] liner permittivities in 
the EPNZ frequency region drive the waveguide even more 
deeply into cutoff. These results affirm that the metamaterial 
liner only significantly alters the propagation characteristics of 
the waveguide in the ENNZ frequency region. 

The complex electric-field magnitudes, obtained using nu- 
merical implementation of the theory developed in Sec. [TTJ 
are shown for both the lined waveguide's forward- wave band 



(Fig. 3(b) ) and its backward- wave band (Fig. 3(c) ), at their 
respective cutoff frequencies. As expected, the forward- wave 
//En -mode cutoff exhibits fields very similar to those at the 
T^ii -mode cutoff of an unlined waveguide. However, in the 
ENNZ regime, two major differences in the field profiles are 
observed: a decrease in field curvature in the inner vacuum 
region and a discontinuous higher concentration of fields 
in the metamaterial-liner region. The latter property can be 
explained using the normal electric-field boundary condition 
which requires that the normal component in the liner must 
increase as |e r 2| — » 0. It should be stressed that the field 
strengths in the vacuum region, although small in comparison 
to those in the liner, are not negligible. Furthermore, the 
reduction in frequency suggests a larger operating wavelength 
inside the vacuum region and, therefore, less field variation. 
It is interesting to note that as e r2 — » 0~ the electric-field 
vectors are collimated, resembling those inside a parallel-plate 
waveguide. In fact, the fields in the liner assume the curvature 
needed to ensure that the tangential component of the electric 
field is zero at the outer PEC boundary. 

IV. Full-Wave Simulations 

Even though the cutoff frequency of the HEu mode can be 
reduced below the unlined waveguide's natural cutoff, it hasn't 
yet been shown if such structures are efficient at transporting 
power. To do so, a transmission analysis is performed using 
the full-wave simulation software HFSS |19] on the simula- 
tion model shown in Fig. |4(a)| In this representative setup, 
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Fig. 4. (a) Full-wave simulation model employed in the transmission 
analysis. A smaller below-cutoff waveguide is placed between two larger 
above-cutoff waveguides and a metamaterial liner possessing the complex 
dispersive permittivity in Sec. [in] is introduced into the smaller waveguide, 
(b) Insertion loss for the unlined case (dashed grey line), the lined case with no 
loss (solid black line) and the lined case with loss (solid grey line). The inset 
shows in detail several resonances in the frequency-reduced backward-wave 
passband. 



two larger vacuum-filled circular waveguides with a radius 
harge — 30mm and a TEu cutoff frequency of 2.928GHz are 
connected by a smaller vacuum-filled waveguide with a radius 
b small — 15mm and a TE\\ cutoff frequency of 5.857GHz. A 
waveport located at the end of one of the larger waveguide 
sections excites the TEu mode at frequencies that lie in 
the propagating region of the larger waveguide, but which 
correspond to the natural evanescent region of the smaller 



waveguide. The dashed curve in Fig. 4(b) presents the insertion 
loss obtained for this setup and verifies that the intermediate 
waveguide under cutoff strongly attenuates the TEu mode. 

Now, a metamaterial liner of thickness t = 1mm is intro- 
duced into the smaller waveguide and assigned the dispersive 
permittivity, e r 2(cj), reported in Sec. Ill According to the 



theoretical results of the previous section, transmission should 
occur below 3. 381 GHz for the frequency -reduced backward- 
wave band and above 5.958GHz in the upper band. Inside the 
frequency-reduced band, the metamaterial liner would effec- 
tively enable a cross- sectional-area reduction of the unlined 
circular waveguide by a minimum of 75%. A 1mm- thick ring 
with the same material properties as the liner is placed on the 
PEC wall at the transition between the two waveguide sections 
in order to aid in the coupling of modes between the differently 
sized waveguides (20). To understand the impact of losses in 
the metamaterial liner on the ability of the lined waveguide to 
transport power, both lossy (u t = 5MHz) and lossless (u t = 0) 
cases are compared. Figure |4(b)| shows the insertion loss in 
the lossless (solid black line) and lossy (solid grey line) cases. 
As predicted by the theory, the metamaterial liner results in 




(c) 



(d) 



(e) 



(f) 



Fig. 5. Complex electric-field magnitudes in the waveguides' H-plane at 
the following transmission peaks in the frequency-reduced backward-wave 
passband: (a) / = 3.381GHz, (b) / = 3.371GHz, (c) / = 3.351GHz, (d) 
/ = 3.324GHz, (e) / = 3.291GHz, and (f) / = 3.256GHz. Each frequency 
corresponds to a resonant condition of an integer number of half- wavelengths 
supported by the lined waveguide section over its length. 



a new passband below / = 3. 381 GHz, corresponding to the 



backward- wave cutoff (/i) in Fig. |3(a)| An upper forward- 
wave band is also observed; however, its cutoff frequency 
is more difficult to infer from the insertion loss due to the 
excitation of high-order modes in the larger waveguide. By 
virtue of the particular geometrical parameters chosen in this 
arrangement (enabling the goal of a 75% cross- sectional- 
area reduction), the HEu mode in the smaller waveguide 
couples to the TMu mode in the larger output waveguide for 
frequencies above 6.09GHz. Another simulation with a smaller 
output waveguide of b\ arQe = 20mm and associated TMu 
mode cutoff frequency of 9.15GHz (not shown) exhibited the 
expected smooth high-pass response. 



The inset of Fig. |4(b)| presents the transmission features 
of the backward- wave region in greater detail, and reveals 
multiple narrow transmission peaks. In both the lossless and 
lossy case, a very fine frequency resolution is required to 
sample the maximum of each peak. When the liner is lossless, 
these peaks achieve total transmission of power through the 
lined waveguide section. Since the fields are concentrated 
strongly in the liner region, the introduction of loss de- 
grades the transmission through the structure; nevertheless, 
the backward-wave band introduced by the liner still exhibits 
a dramatic increase in transmission over the unlined case, 
at times showing enhancements of over 56dB. In fact, even 
in the lossy case, the transmission attains a peak value of 
-21.7dB at / = 3.253GHz. Also evident in Fig. [4(b)] is a 
distinct antiresonance at / = 3.58GHz which is located near 
the plasma frequency of the liner and the edge of the stop 



band (i.e., fa in Fig. |3(a)| ), which appears not to be affected by 
the losses in the liner. The complex electric-field magnitude is 
plotted in the H-plane of the waveguide sections for the first six 



peaks in Fig. 5(a) - 5(f) The E-planes (not shown) exhibit very 
similar field distributions, except that the liner fields are more 
pronounced. Each successive peak corresponds to a resonant 
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(a) (b) 



Fig. 6. Simulated normalized electric field vectors and complex electric field 
magnitudes in the transverse cross-section at the center of the lined waveguide 
at (a) / = 5.958GHz (b) / = 3.38 1G Hz. C lose agreement is observed with 
the theoretical results of Figs. |3(b j| and |3(c)] 




3/k 

Fig. 8. The dispersion curve of the frequency-reduced backward- wave 
passband obtained using three different methods: the full dispersion in 
equation (7a) (solid black line), full- wave eigenmode simulations (solid grey 
circles), and full-wave transmission simulations (empty black circles). 




Fig. 7. Transverse distribution of the Poynting- vector magnitude and direction 
at the center of the lined waveguide section at 3.381GHz. 



condition of an integer number of half- wavelengths supported 
by the lined waveguide section over its length. Moreover, that 
the order of the resonances increases as frequency decreases is 
characteristic of backward- wave propagation, and corroborates 



the dispersion data in Fig. |3(a)| In all cases, it is evident that 
the incoming wavefront is restored at the outgoing side of the 
lined waveguide. These results validate that hollow waveguides 
lined using thin, ENNZ metamaterial liners support unusual 
resonant tunneling phenomena akin to the supercoupling ob- 
served in homogeneously ENZ-filled waveguides. 



Figure |6(a)| presents the transverse electric field vectors in 
the lined waveguide's cross section at the cutoff frequency 
of the upper forward-wave band (i.e., f = fs = 5.958GHz). 
Superposed on these are the complex electric-field magnitudes 
which reveals a T^n-like field distribution that matches well 



with that of Fig. |3(b)| At the backward- wave band cutoff 
of / = 3.381GHz where the first transmission peak is 



observed, Fig. 6(b) verifies the results obtained from the theory 



(Fig. 3(c)): that the transverse fields are strongly confined to 
the metamaterial liner and collimated in the vacuum region. 

It is well known that the direction of power flow can 
vary between different regions in inhomogeneous structures. 
For example, stacked DNG and DPS layers in a parallel- 
plate configuration can transport power forward in the DPS 
region that will, in turn, be coupled backward in the DNG 



region, resulting in zero net power being transported (6j. 
An imbalance of power in each region can lead to non-zero 
net power being directed either forward or backward, as in 
the case of the conventional dielectric-lined waveguides (TJ. 
Figure [7] presents the simulated Poynting vector directions 
and magnitudes at the first resonant peak, along a diameter 
located at the center of the metamaterial-lined waveguide. 
The inset magnifies this behaviour at the interface between 
the vacuum and liner regions. It is evident that power is 
being transported in the positive direction through the liner 
and in the negative direction in the vacuum region at this 
particular location; however, it is also found (not shown) that 
the relative amounts of power in the liner and vacuum regions 
vary continuously along the guide. The total transmission of 
power at these resonant peaks, therefore, suggests a continuous 
resonant exchange of power between the liner and vacuum 
regions. It should be noted that this resonant transmission 
is achieved even though the liner occupies less than 13% of 
the transverse cross-sectional surface area. This allows for the 
waveguide to be miniaturized while still leaving the majority 
of its volume empty for applications which require access to 
its interior. 

To further validate the dispersion of the backward-wave 
band, Fig. [8] compares the dispersion diagrams for the re- 
duced HE n mode obtained using three different methods. The 
first method is based on the theoretically derived dispersion 



relation (7a) and is borrowed from Fig. |3(a)| (solid black 
line). The second method utilizes HFSS's eigenmode solver, 
but employs special techniques to overcome the following 
inherent limitations in the eigenmode solution process: first, 
only propagating modes (i.e. 7 « jf3) may be found; therefore, 
only the backward-wave and forward-wave (and not the com- 
plex) modes can be determined using this method. Second, 
the eigenmode solver determines the resonant frequencies 
of a particular geometry with known boundary conditions 
when the material properties at the solution frequency are 
known a priori (e.g., when the material properties are as- 
sumed constant with frequency). Consequently, in the case 
of dispersive materials, a conventional parametric sweep of 
phase shifts across the structure to obtain the dispersion curve 
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is not a valid approach. However, knowing the frequency at 
which the dispersive metamaterial liner achieves a particular 
permittivity (by way of its known dispersion), the parametric 
sweep may be carried out using this fixed permittivity with the 
knowledge that the obtained dispersion curve is valid only at 
the corresponding frequency. By repeating this process over a 
number of fixed permittivity values corresponding to different 
frequencies, it is possible to obtain the true dispersion curve 
for the metamaterial-lined waveguide from HFSS's eigenmode 
solver. This is indicated in Fig. [8] by the solid grey circles. 
The third method employs the full-wave transmission simu- 
lation results by examining the phase shift through the lined 
waveguide section at each resonant frequency (empty black 
circles). This process has been employed by other research 
groups with experimental results and has demonstrated rela- 
tively high precision, since all the peaks correspond to integer 
numbers of half-wavelength phase shifts [21]. The dispersion 
curves obtained from full-wave eigenmode simulations and 
theoretical analysis correlate strongly with one another across 
the whole backward-wave band, while the data derived from 
the transmission results deviates at higher values of P/k$. This 
is to be expected, since we are comparing the dispersion of 
a finite-length lined-waveguide section to that of an infinitely 
long lined waveguide. 

V. Conclusion 

This work has presented a rigorous study of the propagation 
characteristics of metamaterial-lined circular waveguides. It 
was shown that the introduction of a thin liner possessing 
complex and dispersive permittivity into the waveguide vol- 
ume results in a spectrum of hybrid modes, which may be 
classified as evanescent, propagating, or complex. Whereas 
the fundamental-mode cutoff is largely unaffected by the 
introduction of the liner, a frequency-reduced backward-wave 
passband appears well below cutoff when the real part of the 
liner permittivity assumes negative and near-zero values, and 
it is shown how the cutoff frequency of this mode may be 
designed using the dimensions and permittivities of the liner 
and waveguide. This mode also exhibits uniform, strongly 
collimated fields, with most of the field variation taking 
place in the liner. Full-wave transmission results through a 
lined waveguide section in the backward-wave region reveal 
multiple resonant transmission peaks corresponding to reso- 
nant phase conditions over the waveguide's length. Analytical 
and full- wave- simulation results of the dispersion properties 
and field profiles exhibit excellent agreement, and validate 
that metamaterial-lined circular waveguides operating in the 
frequency-reduced regime offer the potential for miniaturized 
waveguide components suitable for applications in which the 
waveguide volume must remain largely empty. 
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